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Let us consider the matrix Schrddinger equation 

h * E = E* E , h = -D 2 + V (x), D = ^- (1) 

where Vq(x) is an n x n Hermitian matrix with ^-dependent entries, ^ e = 
(ipiE, ■ ■ ■ i V'ti-e)* is a vector of an n-dimensional linear space, and E is a number 
which plays an essential role in different physical applications. For instance, a 
multichannel quantum system is described by this equations pQ. One of the 
most interesting applications of this equation consists in the possibility to involve 
the supersymmetric quantum mechanics [2] for describing scattering of composite 
particles like atom-atom or nucleon-nucleon collisions (SHHEj- In particular, in this 
way one was able to interpret an ambiguity between shallow and deep potentials of 
the nucleon-nucleon interaction [3]. To get a qualitative result authors of Refs jHHHIS] 
apply supersymmetric transformations successively which require performing a lot of 
additional work. In this way one is able to realize only few transformation steps. We 
believe that the progress in applications of this method is essentially delayed because 
of the absence of a simple possibility to get rid of intermediate Hamiltonians and 
go directly to the final result of a chain of transformations. We notice that such 
a possibility exists for the usual (scalar) Schrddinger equation being given by the 
known Crum-Krein determinant formulas [HJ U\ which made it possible to get a 
number of new interesting applications of supersymmetric quantum mechanics in 
describing the nucleon-nucleon scattering Nevertheless, this problem has been 
tackled by Goncharenko and Veselov [Oj when they realized that Gelfand-Retakh 
quasideterminants [TU] may be used for this purpose. We want to point out that 
although their method gives a solution in principle this is is very complicated and 
difficult for practical realization since it involves a matrix calculus with matrices 
defined over a noncommutative ring and in particular it is necessary to invert such 
matrices. 
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In this paper we prove alternative formulas where only usual determinants are 
involved. They are very similar to the known Crum-Krein determinant formulas 
and can be considered as their straightforward generalizations. 



2. First order transformation 



We follow the definition of the Darboux transformation operator given by 
Goncharenko and Veselov [§| defining it as a first order differential operator with 
matrix-valued coefficients 

L = L (x) + Li(x)D (2) 
intertwining ho and hi 

Lh = h x L (3) 
where h is introduced above and hi is defined by the potential V± 

h 1 = -D 2 + V 1 (x) (4) 

If such un operator is found one can get eigenfunctions of hi, by the simple action 
of L on the eigenfunctions of ho, 

$ E = Lty E = (.<PiE,---,<p n E) t , hi<£> E = E<£>E (5) 

Since the equations (JTJ) and (J3J) are homogenous, without loss of generality we 
can put Li equal to the identity, Li — 1 (we suppose detLi ^ 0). After inserting 
L into the intertwining relation (jSJ) we get a system of equations for L and the 
transformed potential Vi. It is not difficult to find the solution to this system [§]. 
Thus, L is given by 

L = D-F } F = U'U- 1 (6) 

and for the potential Vi one obtains 

Vi = V - 2F' (7) 
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The matrix- valued function U is a solution to the equation 

h U = UA (8) 

where A is a constant matrix. 

The known supersymmetric approach PQ-[5] is based on the factorization of the 
Hamiltonian 

h = L + L + XI , A G R (9) 

where L + is defined with the help of the formal relations: D + = —D, i + = —i, 
(AB) + = B + A + and / is the identity matrix. To compare our method with this 
technique let us consider the superposition of L and its conjugate. After a simple 
algebra one finds 

L+L = h - UKW 1 (10) 

In particular, when A = XI, we recover the factorization (jHJ) giving rise to the 
supersymmetry with A meaning the factorization constant. Similarly, the inverse 
superposition gives 

LL + — hi— UNU v (11) 

If A is a diagonal matrix A = diag(Ai, . . . Xn) then the system of equations (jHJ) 
just takes the form of the Schrddinger equation for the columns Uj = (%i, . . . , %,«)' 
of the matrix U = (Ui, . . . , U n ) 

h Uj = XjUj , j = l,...,n (12) 

This means that if we know solutions of the Schrddinger equation (0) then solutions 
of Eq. (JSJ) are also known for the diagonal form of the eigenvalue matrix A. Therefore 
in what follows we will consider only diagonal A's. 

3. Chains of Darboux transformations 



Now we want to consider chains of transformations defined in the previous section. 
Chains appear naturally if we notice that if sufficiently many matrix solutions to 
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the initial equation are known then any such a solution is transformed into a matrix 
solution of the new equation and, hence, the latter may play the role of the initial 
equation for the next transformation step etc. 

Suppose we know iV matrix solutions of Eq. (jHJ) corresponding to different 
eigenvalue matrices A k ^ Ai 

h U k =U k A k , k = l,...,N (13) 

For the first transformation step we take the function U\ and according to © 
construct the transformation operator 

L 1 <_ = D-U[U^ (14) 

We notice that it can be applied not only on vector- valued functions like ^ e but 
also on matrix- valued like U2, • • • , Un- In this way we get the matrix solutions 
V2 = Li^ ^2, • • • , Vn = Li^qUn of the equation with the potential 

V 1 = V -2F{, F 1 =U[U~ 1 (15) 

Now V2 can be taken as transformation function for the Hamiltonian hi = —D 2 + V\ 
to produce the potential 

V 2 = Vi- 2(V 2 V 2 - 1 )' = V - 1F' 2 , F 2 = F x + V' 2 V 2 X (16) 

and the transformation operator L 2< _i — F) — an d so on, till one gets the 

potential 

V N = V - 2F' N (17) 
with Fn defined recursively 

F N = F N ^ + 2Y^Y N 1 , AT =1,2,... , F = (18) 
and Yn being the matrix- valued solution at (N — l)th step of transformations 

YjV = L(jV-l)<-(AT-2) • • • L2^\Li^qUn = L( N _i)^qUn (19) 

which produces the final transformation operator Ljv4-(at-i) — — D + Y^Y^ . 
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To get in this way the final potential Vn one has to calculate all intermediate 
transformation functions performing a huge amount of the numerical work even for 
the scalar case. In practical calculations one is able to perform only few steps which 
restricts considerably possible applications of the method. Fortunately, for the scalar 
case there exists what that are called Crum [H] or Crum-Krein jZj formulas which 
allow to omit all intermediate steps and go directly from ho to h^. The function 

\W(u u . . . ,u N ,ifj E )\ , . 

^ E = ~TfxF( \1 — ( 20 ) 

\W{ui, . . .,u N )\ 

is an eigenfunction of the Hamiltonian with the potential 

provided all Uk, k — 1, . . . , N and ipE are eigenfunctions of the initial Hamiltonian 
ho with the scalar potential Vo: ho = —D 2 + Vo, h Uk = akUk, h ipE = EtpE- Here 
and in what follows the symbol | • | means the usual determinant, W(ui, . . . , Ujv) is 
the Wronsky matrix 

/ Ui U2 ... Un \ 



W(ux, . . .,u N ) 



a, ... u 



N 



(22) 



and the matrix W(ux,...,un) is obtained from W(ui, . . . , Un) by replacing 
its last row composed of u k N ~^ with u k N \ k = 1,...,N. Of course the 
determinant \W(ui, . . . ,un)\ is nothing but the derivative of the determinant of 
the Wronsky matrix |W(iti, . . . , ttjv)| but we write the second logarithmic derivative 
of the Wronskian \W(ui, . . . ,u^)\ in (|2*T|) as the first derivative of the ratio of 
corresponding determinants to stress the similarity of this scalar formula and its 
matrix generalization below. The formula ()20|) really defines for the scalar case the 
superposition of the operators of the type (j!4j) with the replacement of matrix- valued 
functions Z-4 by usual functions 

Ln^o = £jv«-(iv-i) • • • L 2 ^iLi^ (23) 
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ip E = L N ^ Q ip E 



(24) 



We shall prove below generalizations of the formulas (|2Uj) and (121)) to the matrix 
case meaning that we shall solve the recursion defined by (fTHj) . (fTT?|) and find the 
superposition of the operators (|2"3j) but first we need to introduce some new notations 
and to prove an additional statement. 

3.1. Notations 



Define first the nN- dimensional square Wronsky matrix 



W{U 



1) 



Mn) 



u 2 



Z4 



u N \ 



IN 



Here Uk are n x n matrices 



( «l,2;fc 



(N-l) y(N-l) 



U 



(N-l) 
N 



k = l,...,N 



(25) 



(26) 



with columns being n-dimensional vectors Uj-k = (ui t j-k, . . . , w n j ; fc)* so that Uk = 
(Ui-k, ■ ■ ■ , U n -k), k = 1,...,N. We can also present Uk as a collection of rows 
u k = ( u j,i-,k, ■ ■ -Mj> ; fc), 3 = l,---,n,U k = (Ul, Uffi. It will be convenient 
to present (|25j) also in the form 



/ 



W(Ui, . . . Mn-i) 



u N 

Ui 



\ 



N 



y u x ... u N _ x U N 



(27) 



stressing its recursion nature. 
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Introduce also the matrix 

/ 



N) 



W{U u ...,Un) 



(N-l) 



V 



j\(N) 



:28i 



recalling that U° k to be the jth row of the matrix Uk- 
We shall also need the following matrix 

/ 



WHU u ...,Un-i) 



W{U 



i- 



Mn-i) 



Ui-N 

U'. 



\ 



N 



u 



(N-2) 
i:N 



V 



i\(N-i) 



N-l 



.(AT-l) 



U 



(N-l) 



(29) 



First we notice that this is the previous matrix where N is replaced with N — 1 
and ^>e with the vector U^. Another useful remark is that the determinants 
\Wj(Ui, . . . ,Wjv-i)|, i, j = 1, ■ ■ ■ ,n are nothing but minors embordering the block 
W(Ui, . . . , Un-i) in the determinant of the matrix (J27|) . (For the definition of 
embordering minors see Appendix.) 

Finally we introduce the matrices WijiUi, . . . ,Mn), i,j = 1, ■ ■ ■ ,n, constructed 
from the Wronsky matrix (|23|) where the last matrix row composed of matrices 



(N-l) 



is replaced by U k , k 



1,...,N: 



W id (U u ...,U N ) 



u> 



u 



N 



\ 



(N-2) u (N-2) 



u 



(N-2) 
N 

N 



(30) 



The matrices UV 



k> l >3 



1, . . . , n are constructed from the matrix U k 



(N-l) 



by replacing 



its jth row with the it\i row of the matrix IA 



(TV) 



Chains of Darboux transformations for the matrix Schrddinger equation 
3.2. Main Lemma 



9 



In this subsection we prove a lemma we are using in the proof of theorems below. 
Moreover, in proving it as well as the theorems we are using the Sylvester identity 
[TT] which is formulated in the Appendix. 

Consider the matrix 



.4 



a i,p a i,p+i 



\ 



a Pt i . 


(Xpp 




Qp^p+n 




bi,i ■ 


■ hp 


hp+i ■ 


■ &l,p+n 




K hi ■ 


■ hp 


^2,p+l • 


• &2,p+n 


J 



(31) 



Let a be the p x p submatrix of A with the entries aij, i,j = 1, . . . ,p. Denote rrijk 
the minor of A embordering a with jth (j = 1,2) row composed of 6^, j = 1,2, 
i — 1, . . . ,p and kth column (p < k < p + n). Let also m*| be the minor obtained 
from rrijk by replacing its sth row composed of a s j, (s < p) with {p + t)th row 
composed of b^i (t = 1,2). Let now a ts be obtained from a with the help of the same 
replacement, i.e. with the replacement of its sth row composed of a s j, (s < p) by 
(p + t)th row of A composed of b t j (t = 1,2). 



Lemma 1 If \a\ ^ then the following determinant identity takes place 

\ a \ m % = \ ats \ m jk — \oP s \mtk 
Proof. Consider an auxiliary square matrix 



(32) 



/ 



A jt 



Oil 



a>xj> a x,k 



a s ,x ■ 


&s,p 




a p ,x ■ 


^p,p 


&p,k 


b J,i ■ 


■ bj tP 


bj,k 


\ hi ■ 


■ b t , P 





(33) 



Chains of Darboux transformations for the matrix Schrddinger equation 10 

where j,t = 1,2 and the last column containing only two nonzero entries. Take its 
main minor \a\. There are only four minors of Ajt embordering \a\. Minors rrijk 



and \a\ emborder it by the row b 



1, . . . ,p and the next to last and the last 



columns respectively and minors m t k and \a\ emborder it with the last row and the 
same columns. According to the Sylvester identity the determinant composed of 
these embordering minors is equal 



m jk \a\ - m tk \a\ 



\a\\A jt \ 



(34) 



where we can cancel \a\ since it is supposed to be different from zero. 
Now interchange in the matrix Ajt the sth and the last rows to get 



/ 



.4 



jt 



dip a\k 



\ 



hi ■ 


■ b tp 


btk 


1 


dpi . 


(Xpp 







b 3l . 


■ bjp 


bjk 


1 


\ a sl . 


(l sp 




/ 



sth row 



(35) 



The upper-left block of this matrix of dimension p x p is the above introduced matrix 
a ts . Let us find embordering minors for this submatrix. It is clear that m*-^ and 
—m t k emborder it with the next to last column and the the next to last and the 
last rows respectively. The minor embordering a ts with the next to last row and the 
last column has in the last column only two nonzero entries which are equal to one. 
Therefore we can decompose it on this column and after corresponding interchange 
of the rows one gets for it the expression \a ts \ — \oP s \. The last embordering minor 
is equal to — \a\ which becomes evident after corresponding interchange of the rows. 
Once again we consider the determinant composed of these minors and calculate it 
using the Sylvester identity 



\A jt \ \a ts \ 



.ts\ 



As 



m tk \a~-\ -m tk \a JS \ - \a\mf k (36) 
Since \Aj t \ = —\Aj t \ the lemma follows from the equations (f3"4"j) and (p^Hj) . □ 
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In this Subsection we formulate and prove the theorem about the transformation 
of a vector ^e by a chain of transformations introduced at the beginning of this 
Section. 

Theorem 1 The resulting action of a chain of Darboux transformations applied to 
a vector 

®e = (ipiE, • • • , ^ne) 1 (37) 

is the vector 

&e = £j\r<-(jv-i) • • • L 2 ^iL 1 ^ ^ E = (ipiE, ■ ■ ■ , <Pne) (38) 

with the entries ifj E given by 

\W jE (U lt . . . ,U N )\ 

- \w{u^...Mn)\ m 

where WjEipfi, ■ ■ ■ Mn) is given in P?%J) and W(Ui, . . . ,Un) is defined by (j^jl . 

Proof. To prove Theorem 1 we use the perfect induction method. So, first we shall 
prove it for N = 1. In this case 

$ E = L^e = (D- U[U^ E = (<p 1E , • • • , ipnsY (40) 

Denote Ay the cofactor of the element Uj^x in the matrix U\. Then according to 
the definition of the inverse matrix one has 

' ' i=l 

which for the elements of the vector &e implies 

1 n A 
Vie = dipiE -TrnY, u'u^A^e = Tzfr (42) 



*J=1 



Chains of Darboux transformations for the matrix Schrddinger equation 

Consider now the matrix 



12 



IE 



\ 



nE 



V 



u 



u 



j,n;l 



(43) 



If we decompose the determinant on the elements of the last row and all 

determinants appearing in this decomposition except for the one coinciding with \U\ \ 
decompose on the elements of the last column, the resulting expression will coincide 
exactly with the numerator of the right-hand side of (}4*2|) meaning that 

\W jE {U x )\ 



(44) 



which proves the assertion for N = 1. 

Suppose Theorem 1 to be true for the chain of N — 1 transformations meaning 
the following vector gives the resulting action of this chain: 



@£ = -^(Af-l)^O^B = (&1E, ■ ■ ■ , OnEf 



where 



Now since L 



\W jE (U u ...,U N - 1 ) 
\W(U U 



Mn-i, 



(45) 



(46) 



L 



JV-t-(iV-l)MAr-l) 



L(jv-d«-o we have to apply the first order operator 



■/N*-(N-1) 



to the vector 



but first we need to act with L 



(AT-1)< 



_o on the vectors 



Ui-N, i = l,...,n which form the columns of the matrix-valued transformation 



function Um to find the transformation function Y/y = L 
transformation step and determine the operator Ljv<-(at-i) 
supposition of (|4*Hj) . (jjH)) this result is given by 

L(N-i)^-oUi;N = (yi,i, ■ ■ ■ , y n ,if 

\w;(u u ...Mn-i)\ 



(n~i)>^qUn for the A^th 
-D + Y^Y N \ By 



Vj, 



hJ 



n 



(47) 
(48) 



\W(U 1 ,...Mn-i)\ 

This means that the matrix has yj. t (I48j) as entries. If we notice that they 
coincide up to the constant factor 1/\W(JA\, . . . , Un-i)\ with the minors embordering 



Chains of Darboux transformations for the matrix Schrddinger equation 13 

the block W(Ui, . . . Mn-i) in the matrix (J2T|) we can apply the Sylvester identity 
to calculate the determinant |Yjy| which gives 

\W{U u ...,Un)\ 



N\ 



(49) 



\W(U h ...,UN-l)\ 

We need this determinant since the action of the first order Darboux transformation 
operator on a vector is given by (jlUjl . and (f4*Hj) which implies that the vector 



$ E = [D- Y' N Y N l ]<d E 



has the entries 



<PjE 



\Y ne \/\Yn\ 



with 



/ 



1 NE 



l,...,n 



He 



(50) 



(51) 



Yr 



N 







TlE 



\ Vj,i ■■■ Vj,n ".a: ) 
To calculate the entries of this matrix we need to differentiate yj t . 



(52) 

which yields 

\W(U u ...,Un-i)\^ 1 ' \W{U u ...Mn-x)\ (53) 
Here we first calculate the derivative of the determinant \Wj\ keeping in this 
expression the derivative of its last row as a separate term 



\W(U 1 ,...,Un-i)\' \W}{U 1: . . . Mn-x)\' 
•Vj,i + 



\WUU U ...,U, 



N-l, 



(54) 



Here is the determinant of the same matrix (|29jl where only the last row is 
differentiated and A"J is the determinant of the same matrix where in the last 
matrix row of the block W{U\, . . . ,Un-i) one has to differentiate only the mth row 



of the matrices 2 ' ) , k 



N. 



We notice that the matrix of the determinant A™ can also be obtained by 
erchanging one of the rows of matrices uj, N 2 ^ with a row of matrices uj. N 1 \ 
= 1, . . . , N in the minor \Wj(Ui, . . . ,Un-i)\ when it is considered as a minor 
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embordering the block W(U\, . . . ,Un-i) i n the matrix (J27)) . This matrix contains 
two rows with the derivatives of (N — l)th order of elements of matrices Uk in 
contrast with all minors embordering the block W(Ui, . . . ,Un-i) in the matrix (|27p 
which contain only one such a derivative. Therefore one is not able to apply the 
Sylvester identity to calculate this determinant. Just for this purpose we proved our 
main Lemma which implies 

\W(U!, 



= \w mm {u u ...MN-i)\\w}{u 1 , 

-\W jm (U 1 ,...,U N -i)\\Wi(U 1 , 
Now from (JH3|) and (JHHj) one gets 
\Wj(U 1 ,..14 N - 1 )\' 



Mn-i, 
Mn-x)\ 



(55) 



\w(u u ...m. 



N-l, 



{r ji \w{u 1 ,...M. 



N-li 



+ [\W mm (Ui,...,UN-i)\\Wj(U 1 ,...,U 

-\W jm {U x , . . .Mn-i)\ WliUu . . . Mn-i)\]} 
Inserting this into (|53p and taking into account the relation 

n 

\W(U U . . . MN-l)\ = \ W mm(Ul, ---M 



N-l, 



(56) 



W-l, 



m=l 



which is a direct consequence of the structure of the matrices W(Ui, 
W mm (Ui, . . . Mn-i) we get 

1 



N-l 



3,i 



\W(U 1 ,...,U N - l ) 



r ji - y %2\w j , m (Ui,...,u 



N-ljWmA 



m=l 



(57) 
and 

(58) 



By the same means for the vector 0£ (|43j) we can get a similar relation. 

Thus, the determinant of the matrix (|52|) can be written as a sum of two other 
determinants one of which has the last row to be a linear combination of other 
rows and, hence, this determinant vanishes. The matrix of another determinant 
consists of minors embordering block W(Ui, . . . ,Un-i) i n the determinant ((2S1) (up 
to a common factor). Applying once again the Sylvester identity one finally gets 
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which together with (|5U|h (JoTj) and (jjU]) proves the theorem. 
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□ 



3.4- Transformation of potential 

According to (fTTj) to find the potential resulting from a chain of N Darboux 
transformations we have to resolve the recursion defined in (fTHj) and (fTTJj) . This 
is done by the following 

Theorem 2 Let the matrix be defined by the recursion relation F^ = -F/v-i + 

^jv^jv > N = 1,2, . . . with the initial condition F = and Y/v — £(jv-i)<-o^6v with 

the operator L^-i)i~o = L(at_i)^(at_ 2 ) • . . . • L 2 ^\ ■ Li<-o defined in the Theorem 1. 

Then the elements fK of the matrix F^ are expressed in terms of transformation 

functions Uh, k = 1, . . . , iV as follows 

N IW^, . . . ,U N )\ 
J ^ \W{U U ...,U N )\ 1 ; 

where W(U\, . . . Mn) is defined in (p?5j) and Wy(Ui, . . . Mn) is given by (jffi^l . 

Proof. This theorem is also proved by the perfect induction method. For N — 1 one 
has Y\ = 1A\ and W(U\) = IA\. Eq. (ffiU|) follows from inverting the corresponding 
matrix. 

Suppose the Theorem to hold for (N — l)th transformation steps meaning that 
the matrix F^-i has the entries 

f N-i = \W ij {U 1 ,...,U N ^ 1 )\ 

J ^ \W{U u ...Mn-x)\ 1 ) 

To prove the statement we have to calculate the value Fn = Y^Y^ 1 . The matrix 
Ypf has the derivatives (|58J) as the entries. Now since 

(Y N \j = T^Ajt (62) 
\ y n\ 
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where Aij is the cofactor of the element (ijv)i,j in the matrix Y^ then for the entries 
of the matrix Fn one has 

1 1 



f- N 



\W(Ui,...,U N -i)\\Y N 

n 



i=i 



(63) 

FilAji - J] \W im (Ui : . . . ,UN-l)\Vm,lAj 



m=l 



To calculate the first term in the square brackets in (Jfillj) we use the equation 

n 

J2 r u A 3i = \ Y n\ |W(Wi, . . . ,Un-i)\ (64) 



1=1 



where Y$ is the matrix obtained from Yn by replacing in its jth row the 
detrminants \Wj(Ui, . . . ,Un-i)\ (see Eq. I38|) with Tu, I = 1, . . . , n, which follows 
directly from the decomposition of the determinant \Y$\ on its jth row containing 
Fii/\W(Ui, . . . ,ZYjv-i)|- Another relation to be used in these transformations is the 
product of a matrix with its inverse written in terms of matrix elements 

n 

^ VikAjk = 6ij (65) 
k=i 

Now we rewrite equation (|63|) as follows 

7n = \ y n\ _ {Wijiph, ■ ■ ■ ,Un-i)\ / fifi x 
hl \Y N \ \W{U u ...Mn-i)\ 1 ) 

According to (|6*T|) the last term in (JBTIj) represents the elements of the matrix 
Fjv-i- As a result for the entries of the matrix F N one gets 

|yjv| 

The final comment is that the determinants representing numerators of elements 
of the matrix Y^ (we recall that it coincides with the matrix Y^ composed 
of the elements (|48j) except for the jth row composed of the elements 
Fu/\W(Ui, . . . ,ZYjv-i)|) are up to the factor l/\W(Ui, . . . , Wjv-i)| the minors 
embordering the block W(Ui, . . . , Un-i) in the matrix Wy(JAx, . . . Mn) and 
according to Sylvester identity one can write 

^-\W{Ui,...Mn-i)[ { } 

which together with (JBTj) and the expression for \Y N \ proves the Theorem. □ 
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Appendix A 



Here we formulate the Sylvester identity ^1] . Consider a square matrix of dimension 
p + q, p,q= 1,2, ... 



A 



a>ii 



a i,p °i,p+i 



®l,p+q 



\ 



a p,p+q 



ttp t l ■ ■ ■ ^p,p *^p,p+l 

&i,i • • • b\,p &i,z>+i • ■ • ^i,p+9 



(a.i; 



\ bq^l . . . b q ^p &<j,p+l • ■ • ^<J,p+g / 

Let a be the submatrix of dimension p x p composed of the elements ay, i,j = 
1, . . . ,p. If to the bottom of a we add a line of elements fr^i, . . . , bk )P , to the right 
of a we add a column of elements ai iP+ z, . . . , a P}P+ i and the right bottom corner we 
fill with the element bk, p +i, we obtain a square matrix m^\. One says that m^y is 
obtained from A by embordering the block a with kth row and (p + Z)th column. 
The determinant \m,j t i\ is called an embordering minor in the determinant \A\. Since 
k and / can take the values k, I = 1, . . . , q one has q x q embordering minors from 
which one can construct the matrix M = {rrijj). The Sylvester identity relates the 
determinants \M\, \A\ and \a\ as follows: 



|M| 



19-1 



A\ 



(A.2) 
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